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This paper presents a MCO (multi-criteria optimisation) of regional biomass supply chains for the 
conversion of biomass to energy through the simultaneous maximisation of economic performance and 
minimisation of the environmental and social FPs (footprints). The energy supply-chain model contains 
agricultural, pre-processing, processing, and distribution layers. An integrated model, previously 
developed by the authors, for regional biomass energy network optimisation is used as a basis, and now 
extended for simultaneous assessment of the supply-chain performance based on LCA (Life cycle 
assessment). Several total FPs are introduced for “cradle" to “grave” evaluation, which, besides direct, 
comprises also indirect effects caused by products’ substitutions. In the MCO approach, the annual profit 
is maximised against each FP generating different sets of Pareto optimal solutions, one for each FP. With 
this approach the aggregation of different environmental and/or social pressures is thus avoided. The 
results indicate that total FPs enable the obtaining of more realistic solutions, than in those cases when 
only direct FPs are considered. More profitable and less environmentally harmful solutions can be gained 
with significant reduction in total carbon and total energy FPs. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Currently a lot of research is being performed on the topics of 
alternative power/fuel generation [1—3], Examples include fuel 
cells as a technology [4], the utilisation of solar energy [5], wind [6], 
geothermal [7], or tidal energy [8], as well as biomass energy 
applications (bioalcohols [9,10], biodiesel [11], and biohydrogen, 
biogas, bioheat and bioelectricity as well [12-17]). It is important 
that energy systems should be as sustainable as possible [18] and 
based on low carbon thermal processing. 

When considering biomass utilisation, there are important 
issues for ensuring sustainable agriculture, as well as the optimi¬ 
sation of entire bioenergy supply chains [19], Some of these 
important issues are obtaining energy by thermal conversion, 
water, and chemicals’ usage, the economic situation, the world’s 
growing population, climate change, and biodiversity. The distrib¬ 
uted availability of biomass resources, their low density and high 
moisture content, could result in extensive requirements for 
transportation and pre-treatment, and therefore have a significant 
influence on logistics [16,20,21 ]. Another key issue to be considered 
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in relation to exploiting biomass for energy generation is the 
competition with food production. 

Several papers dealing with the optimisation and environ¬ 
mental assessment of biomass energy supply chains have been 
published. For example Zamboni et al. [22] developed a supply- 
chain optimisation model including environmental issues along 
with the economic one, illustrated by a case study of corn-based 
ethanol production system of northern Italy. Mele et al. [23,24] 
developed modelling framework for the design of supply chains 
for sugar and bioethanol production and illustrated it on a case 
study of co-production of sugar and bioethanol from sugar cane in 
Argentina. Gerber et al. [25] presented a methodology to integrate 
LCA (life cycle assessment) in thermo-economic models, illustrated 
by an application of combined synthetic natural gas and electricity 
production from lignocellulosic biomass. 

However, by our knowledge until now no contribution appeared 
that describes general regional biomass and bioenergy supply- 
chain model and incorporates different environmental and/or 
social issues along with the economic one avoiding subjective 
aggregation of different environmental and/or social pressures, and 
besides that includes direct and indirect effects on the environ¬ 
ment. The aim of this contribution was the development of the 
regional supply chains that include important environmental and 
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social FPs (footprints), based on direct and total FPs. For optimisa¬ 
tion, the previously developed mathematical model [16] is used as 
a basis, and has been now upgraded for the simultaneous assess¬ 
ment of FPs. Previous work by authors [16] mainly considered CFP 
(carbon footprint) as an environmental performance indicator. CFP 
stands for a certain amount of gaseous emissions relevant to 
climate change, and is associated with human production and 
consumption activities [26]. Although the synonyms “climate 
footprint” [26] and “GHG (greenhouse gas) footprint” [27] are 
perhaps more appropriate, the term “CFP” is used in this paper, 
mainly due to its broader acceptance so far. “CFP” is usually defined 
based on the amount of CO2 and other GHGs emitted over the full 
life cycle of a process or product, and expressed in mass unit per 
functional unit [28,29], 

Besides GHG emissions, several other negative impacts may 
result from biomass production and use for bioenergy generation — 
water pollution and shortage, as well as food and land scarcity, to 
name a few. 

Therefore, other environmental FPs should also be considered: 

- Energy footprint — EFP. One of its definitions is that it repre¬ 
sents the demand for non-renewable energy resources [30], 
expressed in energy unit per functional unit; 

- Water footprint - WFP, defined as the total volume of direct 
and indirect freshwater used, consumed and/or polluted [31 ], 
expressed in volume unit per functional unit; 

- Agricultural land footprint - LFP, the agricultural land area 
used for growing biomass for both food and energy [32], 
expressed in area unit per functional unit; 

- Water pollution footprint — WPFP, the amount of substances 
emitted to water in the environment, expressed in mass unit 
per functional unit [33]. 

It should be noted, that environmental FPs are usually consid¬ 
ered to be measured in units of area. However the data expressed in 
units of area show high variability and the high possible inaccuracy 
in results since they would have to be based on a variety of different 
assumptions (see e.g., Ref. [26]). Converting some of the FPs to an 
area units can prove to be problematic, especially for processes that 
are not primarily area-based such as chemical processes [34], 
therefore more appropriate unit for each of considered FP is used. 

Food versus fuel competition for biomass utilisation is another 
very important issue relating to the usage of biomass for fuels (first 
generation biofuels), which should be considered. This problem is 
emphasised by the Nestle chief executive: “If, as predicted, we look 
to use biofuels to satisfy 20% of the growing demand for oil products, 
there will be nothing left to eat" [35], Another problem related to the 


biofuel industry is the global increase of food prices as a result of 
using more crops and land for energy purposes [36,37], For this 
reason, a FEFP (food-to-energy footprint) should also be included as 
a special social FP in this study. It is defined as a mass-flow rate of 
food-intended crops converted into energy, expressed in mass unit 
per functional unit. 

On the other hand, when considering the direct effects that 
different FPs have on the environment, this may result in 
misleading solutions. A much broader picture, leading to more 
complete estimates, can be obtained if the indirect effects caused by 
product substitution are considered, too. For instance, although the 
production of biofuels can cause a possibly significant direct burden 
on the environment, this effect is more than balanced out by the 
fact that biofuels substitute equivalent amounts of more harmful 
fossil-based fuels. Therefore, it is important that FPs account for 
both direct and indirect effects. The concept of total FPs is now 
introduced and applied to the biomass and bioenergy supply chain. 

2. Implementation of sustainable criteria through total FPs 

Biomass and bioenergy supply chains involve harvesting, 
storage, pre-treatment, conversion steps, distribution and usage of 
products, and transport amongst their activities. The superstruc¬ 
ture of the considered supply chains is presented in Fig. 1. 

It contains several layers dedicated to key activities: 

(i) “Harvesting and supply" layer (layer 1 — LI) dedicated to raw 
material supplies 

(ii) “Collection and pre-processing” layer (layer 2 - L2) for the 
pre-treatment of biomass when obtaining intermediates 
directly suitable for energy generation 

(iii) “Main processing" layer (layer 3 — L3) comprising the core 
processes for energy and other bio-products’ generation 

(iv) “Use” layer (layer 4 - L4) comprising the usages of the energy 
products. 

The transportation activities have no dedicated layer. Rather, 
they are embedded in the communication arrows between the 
layers. 

The bioenergy supply-chain model, which consists of mass 
balances, production and conversion constraints, cost functions, 
profit objective function, and CFP calculation [16,17] has been 
extended to include the consideration of various environmental 
FPs. In addition to these, a FEFP is defined in order to evaluate the 
risk of diverting farmland for the production of fuel rather than 
food [38], As a result, the mathematical model employs important 
direct, indirect, and total FPs (the sum of direct and indirect FPs) 
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Fig. 1 . The superstructure of the biomass and bioenergy supply chains (after Ref. [16]). 
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along the whole supply chain. Different environmental FPs and 
FEFP are defined below, as per year and per unit of the supply-chain 
network’s total area (71, km 2 ). 

The model is formulated in the MI(N)LP (mixed-integer (non¬ 
linear programming) form. The following sets have been defined 
within the model (Eqs. (1)—(9)): 

i) Set / for supply zones within harvesting and supply layer (LI) 
with elements i e I 

ii) Set M for collection and pre-processing centres within 
collection and pre-processing layer (L2) with elements bieM 

iii) Set N for processing plants within main processing layer (L3) 
with elements n e N 

iv) Set J for demand locations within use layer (L4) (see Fig. 1) 
with elements j e J 

v) Set T for process technology options with elements t e T 

vi) Set P for products with elements p e P 

vii) Set PI for intermediate products PI E P with elements pi e PI 

viii) Set PD for directly used products (PD £ P) from the layer L2 
with elements pd g PD 

ix) Set PP for products coming from the layer L3 (PP £ P) with 
elements pp g PP 

x) Set PS for substituted products with elements ps g PS 

xi) Set F for environmental and social FPs with elements fe F 

xii) Set PT=PIxT for pairs of intermediate product, and the 
applicable process technology for it with elements (pi,t) g PT 

xiii) Set PSP=PS xPP for pairs of substituted product and 
produced product (if a produced product substitutes 
conventional product) with elements (ps.pp) g PSP 

xiv) Set PIP = PIx PP for pairs of intermediate products and 
produced products (if a produced product is produced from 
a given intermediate product) with elements (pi,pp) g PIP 

The definition of a direct environmental FP is: 

• For the supply layer, ENVB|' is calculated from the production 
rate q-" p! LI of biomass type pi g PI in supply zone i e I, multiplied 
by the specific environmental FP eip//,/e F, caused by growing 
the biomass. The result is further summed over all the supply 
zones and biomass types, and divided by the supply network’s 
total area (A): 


- 


or the pre-treatment layer, ENVBp - the mass-flow rate of 
liomass pi from the supply zone i in LI to the pre-treatment 
entre me Min L2, multiplied by the specific environmental 
P for that biomass, caused by pre-processing, and also 
valuated per unit of total area (A) after the summation over all 
ones, pre-treatment centres, and biomass types: 




types, technology options and pairs of intermediate products 
and the applicable process technology for it, divided by the 
network’s total area (A): 


envb- = (e E E E 

\neJVpieP7<=P tel (pi,t)ePT 


For the use layer, ENVB| 4 is calculated as the mass-flow rate of 
the resources directly used from the collection and pre¬ 
processing layer (Fig. 1) pd g PD and also those products 
coming from the processing layer (plants) pp g PP, each 
multiplied by the specific environmental FP for this product, 
eip^ and eip 4 ^, caused by their usage. The summation is per¬ 
formed over all pre-treatment centres (for direct products) or 
processing plants (for converted products), demand locations 
and product types per network’s total area (A): 


EE 

+ £SJk«‘-*«')A v/sF <4> 

Finally, for the transportation of materials between layers, the 
environmental FP ENVBf is also defined per network’s total 
area (A). Its value depends on the density of the biomass, 
distances (D^ Lb ), mode of transport, rate of biomass supply 
($£“). and road conditions (/™ d ’ LaLb ). The specific envi¬ 
ronmental FP is defined per t km: 


ENVB} 4 = 


EEJ c CC d,L1 "' ei ?r’C 

-e&r - 45 S 3 

+JC E pUl D ”^ L4 ■•C d,L2,L4 ' ei pdf L4 ' *® L4 
+EEE s C 4 C dW 4 -C u -C 

/A V/eF 


Eq. (5) represents the environmental FP caused by transporting 
biomass pi from LI to L2 (the first term), intermediate product pi 
from L2 to L3 (the second term), the directly used product pd from 
L2 to L4 (the third term), and the converted product in plant pp 
from L3 to L4 (the last term). 

Note that specific environmental FPs are defined in t/MJ for heat, 
t/MWh for electricity, and t/t for other products. 


For the processing layer, ENVB| 3 is obtained as the mass-flow 
rate of the intermediate product pi to the selected technology 
te Tat process plant n e N in 13, multiplied by the specific 
environmental FP for that product, eij^t/, caused by processing 
and summed over all process plants, intermediate product 


The indirect environmental footprint represents unburdening by 
product’s substitution. When conventional energy is replaced 
by biomass energy, the indirect effect is a reduction of the FP. 
The indirect environmental FP for the use of biomass is then 
defined as the unburdening related to the substitution of 
conventional, mainly non-renewable products. The correlation 



L Cucek et at / Energy 44 (2012) 135-145 


between conventional energy and biomass energy is defined 
by the substitution factor/^*,, where ps e PS is the substituted 
product, and pp e PP is the produced biomass product. With 
the generation of bioenergy and bioproducts, the supply and 
consumption of previously used resources is reduced in the 
amount specified by the substitution factors. The indirect 
environmental FP is defined as the product between the 
specific environmental FP for the substituted product, eip" d , the 
substitution factor for that product, /p“p p , and the mass-flow of 
the biomass product converted in the plant, j m 14 , per unit of 
total area. The summation is performed over all the substituted 
products, produced products and pairs of substituted products, 
and the produced products (if a produced product substitutes 
a conventional product) (first term), and over process plants, 
demand locations, and produced products (second term). 

Note that the indirect environmental FP ENVB| nd is negative, 
since the substitution causes unburdening: 


Firstly, MI(N)LP-1 is used to obtain an initial or reference solution, 
representing the maximum profit possible/upper bound of the 
profit. This is then followed by solving MI(N)LP-2 for obtaining the 
multi-objective optimal solution(s). 

MCO (Multi-criteria optimisation) is performed at the second 
MI(N)LP level (MI(N)LP-2). MCO, also known as multi-objective 
optimisation (MOO), is the process of simultaneously optimising 
two or more in many cases conflicting objectives, subject to certain 
constraints. Different methodologies can be applied for solving 
MOO problems. In this contribution the e-constraint method [39] is 
applied and a sequence of constrained single-objective (MI(N)LP- 
2), problems solved for each FP/se F, as the maximisation of the 
profit P subjected to a relative FP (RFP/). RFP is defined as the FP 
obtained at MI(N)LP-2 (FP(x,y)/) divided by the FP obtained at MI(N) 
LP-1 (FP/). It decreases sequentially from the maximal FP obtained 
at MI(N)LP-1 (RFP— 1) by a suitable step size Ac,/ until there is no 
feasible solution. 

The synthesis problem at MI(N)LP-2 takes the following form: 


ENVB? nd = - 


( E E E ei ps/'/ps U pi 

\pssPS ppePPcp (ps.pp)ePSP 


EE E 

neN jej ppePP<=J 


/A V/cF 


The total environmental footprint of the supply-chain’s network 
is defined as a summation of all the direct FP elements, and the 
indirect FP correction: 


max P = c T y+/(x) - £ ICH-ENVBa. - £ 10- 6 FEFP fr 
** fr*f fr*f 

s.t. Ay + h(x,y) = 0 

By + g(x,y)<0 (MI(N)LP,) ; 

RFP (x,y)f < e t j 

(^° <X<X VP )eX(-R n . y e y={0,l} m 
e i 11 / = e if ~ Ae if 
where RFP is defined as 


ENVBf = ENVBr + ENVBp + ENVBP + ENVB? 4 

■ * (7) 

+ ENVB} r + ENVB} nd 

• The FEFP is only defined for those multi-functional (multi¬ 
product) crops that can result in a supply of food, fodder, and 
energy. Biomass crops, which may take land that could be used 
for growing food crops, are not considered in this FP (FEFP is 
zero). This social FP is defined as a mass-flow rate of energy 
converted from food-intended crops. The FP is summed over all 
process plants, intermediate products going to plants, con¬ 
verted products from food crops in plants, technological 
options, and pairs of food-crop products and the applicable 
process technology for them (numerator), and over all the 
supply zones, pre-treatment centres and biomass types 
(denominator). Weight loss by drying and the conversion of 
intermediates into products is described using the conversion 
factor/p“ nv,L2 . The FEFP is divided by the network’s total area: 


fefp=(e e e e e 

yneNpieP/cpppePPcp t eT (pi,pp)eP/P 

>)))A 

(8) 

2.1. The two-level multi-criteria approach 

The synthesis is performed in two levels. At the first MI(N)LP 
level (Ml(N)LP-l) the synthesis model from the previous work 
[16,17] plus different environmental FPs (ENVB/) and FEFP 
(Eqs. (1)—(8)), is solved with the objective of maximising profit. 


RFP (x,y) f = 


FP(x,y)| 

FP? 


O) 


x denotes the vector of continuous variables (flow-rates, design 
variables, etc.), y denotes the vector of binary variables (discrete 
decisions), and/(x), h(x,y) and g(x,y) are the continuous non-linear 
constraint functions. 

Note that MINLP/is performed separately for each FP/,/e F, while 
all the rest FPs Jr.frc FJr 3= F, are minimised with a small weight to 
provide solutions with the least values for those FPs in cases where 
multiple FP solutions exist. The weight should be very small in order 
not to interfere with the maximisation of profit, e.g., 10 -6 . 

At Ml(N)LP-2 two sets of MCOs were performed: i) maximising 
the profit whilst simultaneously minimising the relative direct FPs 
and ii) maximising the profit whilst simultaneously minimising 
the relative total FPs (direct + indirect FPs). Different sets of Pareto 
optimal solutions, one for each FP, are generated at MI(N)LP-2 as 
2-dimensional projections of a multi-dimensional problem, and 
so the subjective aggregation of different FPs, since they are 
usually expressed in different units, is thus avoided. Note that in 
case of negative FPs e-inequality constraint is changed into an 
equality constraint and a single-parametric optimisation is per¬ 
formed. In this way, non-trade-off optimal solutions are obtained 
by the sequence of (MINLP-2),- where profit is maximised while 
relative FPs are incrementally increased from -1 to zero of even 
some positive values until there is no feasible solution at zero 
profit. 


3. Demonstration case study 

The above concept and developed model was applied within 
a regional biomass and bioenergy supply chain [16] (see Fig. 2). The 
region was divided into 10 zones (supply layers) with a total area of 
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1000 km 2 . In the region there were six optional collection centres 
(m), three optional plants (n), two locations of local (j°) and °ne of 
cross-regional demand (f). The geographical features of this 
network considered in the case study, are illustrated in Fig. 2. 

Several biomass types, technology options and bioproducts 
were considered in the synthesis. The biomass types utilised were 
corn, corn stover, wood chips, MSW (municipal solid waste), 
manure, and timber. The technologies for converting biomass were 
the dry-grind process, anaerobic digestion, MSW incineration, 
timber sawing and incineration, and the produced bioproducts 
were bioelectricity, bioheat, bioethanol, boards, digestate, and 
DDGS (distiller’s dried grains with solubles). Environmental impact 
was assessed through CFP. This case study [16] was now extended 
to employ different FPs originating from bioenergy, except CFP, 
which had already been evaluated in previous work [16,17], 

The two newly-defined sets were: set F for FPs, F — {CFP, 
EFP, WFP, LFP, WPFP, FEFP} and set PS for substituted products, 
PS = {electricity mix, heat mix, gasoline, corn, inorganic fertiliser, 
steel}. 

The data for various specific FPs were taken from many different 
sources, and are given in Tables 1—3 for direct, and in Table 4 for 
indirect FPs. The rest of the data can be found in Cucek et al. [16]. 

3.1. Direct FPs 


The data for raw material’s various direct specific environmental 
FPs ei are given in Table 1. Specific environmental FPs are defined in 
t (WFP, WPFP and CFP) and in GJ (EF) per t of biomass. 

The data for various direct specific FPs caused by pre-treatment 
and processing per ton of intermediate product; for pre-treatment — 
biomass transported to collection and pre-processing centres, and 
for processing - pre-treated transported intermediate products, are 
given in Table 2. 

In Table 3 direct specific environmental FPs for transport are 
presented. Biomass was transported by truck from the supply 
zone to pre-treatment centres (CC), with consumption of 0.3 L of 
petrol/km. Intermediate products were transported by rail from CCs 
to plants. Transport of direct products (corn intended for food) to 
consumers is assumed to be partly by truck, partly by rail, and 
produced products to consumers by truck. Bioethanol was trans¬ 
ported to filling stations by tank trailers (cisterns). Electricity was 
transported by alternating current power lines; heat was used for 
buildings and water heating. 

The environmental FPs for products are only assumed for 
digestate: CFP 0.017 t/t and WPFP 2.01 kg/t. The data in the litera¬ 
ture vary significantly, environmental burdens are namely changed 
regarding the location, time, technology, composition etc., and, 
therefore, proper evaluation of LCA without measuring the FPs is 
impossible. Also, all other parameters vary considerably, such as 


Direct FPs represent direct burden on the environment when 
biomass is used, pre-treated, processed and biomass products are 
utilised. Burdens are caused also by the transporting of raw mate¬ 
rials to collection and pre-processing centres, to plants or directly to 
consumers, and by the transporting biomass products to consumers. 


Table 1 

Direct specific environmental footprints for biomass. 


Corn grains 

Organic manure 
Wood chips' 
MSW 


Water (t/t) Non-renewable Emissions to 

energy water (t/t) 

consumption (GJ/t) 

900 [40,41] 1.73 [43,44] 0.032 [46] 

900 [40,41] 1.73 [43,44] 0.032 [46] 

0.75 a 

2500 [42] 0.75 [45] 

0.229 a 

1500 [42] 0.5 [45] 


Emissions 
to air (t/t) 


0.154 [43] 
0.154 [43] 

0.066 [47] 
0.044 [47] 


Table 2 

Direct specific environmental footprints for pre-treatment and processing. 


Water (t/t) 


Non-renewable 

energy 

consumption (GJ/t) 


Emissions to 
water (t/t) 


Dry-grind 

process 

MSW 


1.3 [48] 2.5 [53] 

0.3l b [16] - 


0.0016 [55] 


Anaerobic 0.091 [50] 

digestion 

Sawing 10.6 [58] 

Corn stover 0.005 [58] 

compressing 
Corn grain 0.5 [58] 

drying 

Timber 0.004 [58] 

drying 


0.036 [57] 
0.01504 [51,59] 

1.251 [52,56] 
0.0108 a [54,60] 


Emissions 
to air (t/t) 

0.147 [16] 
0.415 [16] 

0.00125 [49] 
0.00262 [59] 

0.09 [16] 
0.00078 [16] 


a Since organic manure and MSW are waste, there is no available data on WFP, 
and only raw materials moisture is taken as WFP. 
b Assuming that from the entire plant mass around 50% is grain and 50% stover [9], 
c Assuming the density of wood chips is approximately half density of timber [16], 


a Air and kiln drying of wood is assumed where thermal energy for kiln drying is 

b A little higher water footprint than moisture content of MSW is assumed 
(22.9%, [16]). 
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Table 3 

Direct specific environmental footprints for transport (in t/(tkm) and in GJ/(tkm)). 


Water (t/(tkm)) [54] 


Manure, digestate, DDGS 

Wood chips 

MSW 

Bioethanol 


a Rail: 0.15 MJ/(m 3 km), road: 


1.36 x 10“ 4 
2.33 x 10 -3 


5.6 x 10 4 
1.24 x 10 4 
2.45 x 10 4 


0.28 MJ/(m 3 • km). 


Rail 

7.3 x 10~ 5 
1.25 x 10‘ 3 

2.63 x 10- 4 


1.31 


10 4 


Non-renewable energy 
consumption (GJ/(tkm)) [45] a 



Emissions to 
water (t/(tkm)) 
Road Rail 


Emissions to air (t/(tkm)J 
[16] 

Road Rail 

5.3 x 10~ 5 8 x 10~ 6 

1.1 x 10~ 3 8 x 10~ 6 

5.3 x 10~ 5 

2.4 x 10 4 8 x 10“ 6 

1.3 x 10~ 4 
2.7 x 10“ 5 

5.3 x 10“ 5 8 x 10~ 6 


biomass availability in the region, investment and operating cost, 
prices for raw materials and products, to name a few. It should be 
noted that this case study is just a demonstration case study which 
can be extended to any specific biomass energy supply chain. 
However, the quality of a solution depends on non-convexities 
involved in the model as when solving complex problems global 
optimality of solutions cannot be guaranteed. 

3.2. Indirect FPs 

Indirect FPs are set of FPs that indirectly unburden or benefit the 
environment when biomass products are utilised instead of 
conventional mainly fossil energy since its direct harmful impact on 
the environment is thus avoided. The correlation between 
conventional energy and biomass energy is defined by substitution 
factor/p S U p p . Substitution factors of substituted products to biomass 
products and their specific environmental FPs are presented in 
Table 4. 

4. Results and discussion 

Applying data from Tables 1^4 in the MI(N)LP supply-chain 
model [16,17], optimal solutions were obtained. At the M1NLP-1 the 
optimal solution of 34 M€/y was obtained and FPs were calculated 
from the optimal solution. 

The integrated MINLP model of regional biomass supply chain 
consisted of around 1100 continuous variables, 20 binary variables 
and 650 constraints and was solved in a fraction of a second. The 
MINLPs were performed using D1COPT/GAMS on a computer with 
2.33 GHz Intel® Core™2 Quad Q8200 processor with 4.00 GB 
of RAM. 

The calculated direct, indirect and total FPs at MINLP-1 are 
presented in Table 5. 

Those values obtained were set as maximal FPs, and when 
normalised, set to 1. 

From Table 5 it can be seen that by replacement of biomass 
energy and bioproducts for conventional energy and products, 
CFP and EFP are reduced, while on the other hand WFP and WPFP 
are increased. LFP and FEFP are zero for conventional fossil 


energy and therefore total LFP and FEFP remains the same as for 
direct FPs. 

At MINLP-2 simultaneous maximisation of the profit and min¬ 
imisation of the i) relative direct FPs and ii) relative total FPs were 
performed. By maximising the profit at MINLP-2, for each FP/out of 
six of them, a sequence of (MINLP,)j was performed where a relative 
FP RFP/was decreased from 1 (also increased from -1 for relative 
total FPs) to 0 (some positive value for relative total FPs) whilst the 
other FPs ff, were calculated. Thus a set of Pareto curves was 
obtained, one for each FP [71], 

4.1. Direct FPs 

In this paper only the main Pareto curves are presented where 
a normalised FP/was decreased from 1 to 0 by suitable step size. In 
Fig. 3 profit versus relative direct FPs are presented. Note that FPs 
are normalised so that each FP is divided by its maximal value 
obtained at MINLP-1 (Eq. (9)). 

It can be seen that the FEFP especially differed from FPs. The 
lower profit was obtained by lower FEFPs, and that means that it is 
economically more preferable to produce energy rather than food. 
All other FPs showed similar behaviour with slight concave 
curvatures caused by the fact that for increasing the profit less and 
less sustainable alternatives should be selected. The profit 
increased at the lowest RFPs most steeper, for approximately 
4.7 M€/y per 0.1 of RFP, while at the RFP 0.3 and higher only for 
approximately 2.9 M€/y per 0.1 of RFP. 

Heat and electricity were generated by all ranges of different 
FPs; bioethanol and DDGS as well were produced by all ranges of 
RFPs except by lower CFPs (below 0.25) and EFPs (below 0.65), 
anaerobic digestion was selected within the whole range of EFPs, 
LFPs, WFPs, and FEFPs, and in the lower range of CFPs and WPFPs. 
Board-making was selected by the whole range of LFPs, WPFPs, 
and FEFPs, and food was not produced by the whole range of 
WFPs, WPFPs and LFPs. In Fig. 4 profit versus relative products 
flow-rates are presented when relative direct CFP was simulta¬ 
neously minimised with maximising the profit. The maximal 
flow-rates of products when their relative flow-rates are 1, are: 
243 GWh/y of electricity, 1.25 PJ/y of heat, 62.2 kt/y of ethanol, 


m factors and specific environmental footprints of subsl 


Substitution factor 


tctricity mix/electricity 


1/1 


Heat mix/heat 
Gasoline/bioethanol 1/1.53 

Com/DDGS 1/1.2 [70] 

Inorganic fertilizer/digestate 1/34 [61,62] 

Steel/Boards 1/1 


0.175 t/GJ [58] a 
3 t/t [67] 

900 t/t [66] 

150 t/t [65] 

260 t/t [65] 


3.6 GJ/MWh 
1 GJ/G] 

44.4 GJ/t 
15.9 GJ/t [16] 
40 GJ/t [44] 
35 GJ/t [64] 


0.12 [68,69] 
0.032 [46] 
1.15 [68,69] b 
0.02 [68,69] 


0.15 t/MWh [63] 
0.076 t/GJ [63] 
3.9 t/t [63] 

0.154 t/t [43] 
2.508 t/t [63] 
1.609 t/t [63] 


a The average WFP of a i 
b The values for the amounts of used fertilisers ar 


x of ur 


:n from Ref. [68] for corn, production in the USA. 
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Table 5 

Direct, indirect and total footprints obtained at MINLP-1. 



Direct footprints 

Indirect footprints 

Total footprints 

CFP (t/(km 2 y)) 

117.65 

-311.95 

-194.3 

WFP (t/(km 2 y)) 

376,500.75 

-39,210.75 

337,290 

EFP (GJ/(km 2 y)) 

1440.65 

-4906.72 

-3466.07 

WPFP (t/(km 2 y)) 

12.02 

-6.47 

5.55 

LFP (km 2 /(km 2 y)) 

0.32 

0 

0.32 

FEFP(-) 

0.38 

0 

0.38 


3.65 kt/y of organic fertiliser, 48.1 kt/y of ethanol, 9.5 kt/y of 
boards and 192 kt/y of com. 

4.2. Total FPs 

Total FPs are composed of direct and indirect FPs. The values at 
MINLP-1 (last column in Table 5) were used to normalise FPs 
obtained at M1NLP-2. Total FPs were divided by their absolute value 
of FP at MINLP-1 and therefore relative total FPs were presented on 
the x-coordinate from +1 to 0 (positive total FPs) and from -1 to 
0 or even higher value (negative total FPs). In Fig. 5 Pareto curves for 
profit versus relative total FPs are presented. In case of initially 
negative FPs the curves are obtained by the sequence of (MINLP-2); 
where profit is maximised while relative total FPs are incrementally 
increased from -1 to some positive values until there is no feasible 
solution at zero profit. In those cases the MCO changed to single- 
parametric (MINLP-2)j, where relative FPs RFP/are equal to e, and 
sets of non-trade-off optimal solutions are thus obtained, since 
relative FPs lowered despite that profit increased. In case of initially 
positive FPs the profit is maximised while their relative total FP is 
decreased incrementally from +1 to 0. 

From the above Fig. 5 it can be seen that total CFPs and EFPs 
showed similar behaviour (mainly negative), from which it can be 
gathered that they are in relation. Similar behaviour also showed 
WFPs, WPFPs and LFPs and are always positive. The reasons are: 
biomass energy requires more water, large land areas and some 
biomass also requires large amounts of chemicals which cause 
pollution to water. All this is not the case for conventional energy 
(mainly fossil). From Fig. 5 it can be seen that with biomass energy 
CFPs and EFPs are reduced compared to currently used energy, 
while on the other hand WFPs, WPFPs, LFPs and FEFPs are 
increased. It also should be noted that by zero total CFP 


Direct footprint < 



— Carbon footprint, CFP 

— Water footprint, WFP 

— Food-to-energy footprint, FEFP 


• Energy footprint, EFP 

-m -Water pollution footprint, WPFP 

-1-Land footprint, LFP 



much higher profit (14.4 M€/y) is obtained compared to direct 
alternative. 

The production of biomass products versus total FPs shows 
similar trends as versus direct FPs. In Fig. 6 relative products flow- 
rates versus profit are presented when total relative CFP was 
simultaneously stepwise increased to zero profit with maximising 
the profit. The maximal flow-rates of products are the same as by 
direct FPs. Heat and electricity were again generated by all ranges of 
different FPs, however their production was higher compared to 
the direct alternative, by the zero CFP 145 GWh/y of electricity and 
0.75 PJ/y of heat were generated (1.2 kWh/y of electricity and 
6.03 MJ/y of heat for direct alternative), and by the zero profit 
35.3 GWh/y of electricity and 0.18 PJ/y of heat were generated. 
Bioethanol, DDGS and boards showed similar trends for relative 
direct and total CFPs. Corn production was increased with 
increasing CFP, by the zero profit maximal amount of corn was 
produced, and digestate was only produced by the lowest CFP. 


Total footprint < e 



-A- Carbon footprint, CFP •- Energy footprint, EFP 

-*- Water footprint, WFP ■- Water pollution footprint, WPFP 

-X- Food-to-energy footprint, FEFP 1-Land footprint, LFP 


Fig. 5. Profit versus total 


Fig. 3. Profit ■ 


; direct specific relative footprints. 


footprints. 
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Total carbon footprint < £ 



-*- DDGS -*- Boards -1-Com 

Fig. 6. Profit versus relative products flow-rates for the case where total CFP was 

stepwise increased to zero. 


4.3. Total/direct FPs 

If relative FPs for total FPs are defined as total FPs divided by 
values of the corresponding direct FPs obtained at MINLP-1 (rela¬ 
tive total/direct FPs), then we get some additional information 
about how much better (or worse) are total FPs than direct FPs 
because of product’s substitution effects. If a relative total/direct FP 
evaluated at maximal value of profit gets a value less than 1 means 
that a total FP has a lower value than its corresponding direct FP 
and consequently, product’s substitution helps decreasing envi¬ 
ronmental burdening, especially when it gets very negative values. 
Fig. 7 shows Pareto curves for profit versus relative total/direct FP 
for the case study. All relative total/direct FPs but for LFP and FEFP 


are decreased compared to their direct impacts, especially EFP and 
CFP which are very negative and substantially unburden the envi¬ 
ronment compared to currently used energy. Positive but less than 
1 are relative total/direct FPs for WPFP and WFP meaning that their 
burdens are decreased compared to their direct impacts but not 
diminished. Only relative total FPs for LFP and FEFP are the same as 
their direct impacts, because production of biomass requires large 
land areas and production of energy from food-intended crops 
consumes large amounts of these crops. Note that absolute values 
of total FPs at maximal profit are written at the ends of their Pareto 
curves. Note also that these total FP curves in Fig. 7 can be now 
directly compared to direct FP ones in Fig. 3 since all of them are 
normalised using the same FP values. 

5. Conclusions 

This contribution presented MCO of the regional biomass and 
bioenergy supply chains in order to maximise the economic 
performance of the supply chain by simultaneously minimising 
different environmental FPs and special FEFP. Direct, indirect 
(compared to the current situation) and also total (direct + indirect) 
environmental FPs were studied in this multi-criteria approach. The 
problem was solved in two steps, at MINLP-1 profit was maximised 
and ‘maximal’ (‘minimal’ total CFPs and EFPs) FPs were calculated, 
and then at MINLP-2 MOO was performed, where simultaneously 
profit was maximised and each RFP was stepwise minimised 
(stepwise increased from -1 for CFPs and EFPs). Different trade-offs 
were obtained between the network’s alternatives. 

From the optimisation of the total FPs it was seen that CFP and EFP 
are reduced while WFP, WPFP and LFP are increased compared to 
conventional mainly fossil energy. Biomass energy compared to 
conventional energy requires much more water, transport, chem¬ 
icals, cause pollution to water, and it requires also large land areas. It 
is interesting to note that CFP and EFP on the one hand, and WFP, 
WPFP and LFP on the other hand showed similar common behaviour. 
It is clear that product substitution can have significant effect on 
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-•- Water pollution footprint, WPFP -X- Food-to-energy footprint, FEFP -1-Land footprint, LFP 

Fig. 7. Profit versus relative total/direct footprint. 
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total FPs and should be considered. According to the solution 
obtained by minimising relative total/direct FPs, total EFP is reduced 
by about factor -2.4, total CFP by -1.4, and even relative total WPFP 
and WFP to +0.5 and +0.8, when compared to their direct values. 

In the future the number of environmental FPs should be 
reduced to a minimum of “independent” FPs through correlations 
among FPs. In the set only one FP, an “independent” one, can be 
considered within the MOO whilst the rest of the “dependent” FPs 
can be evaluated after optimisation from the independent ones. In 
this way, the dimensionality of the criteria within the MOO can be 
significantly reduced, so that a multi-parametric optimisation can 
be performed with independent FPs as parameters, thus avoiding 
the subjective weighting of environmental pressures. The correla¬ 
tions from other FPs should also be investigated, as assessed in this 
contribution, such as nitrogen and phosphorus FPs, and the issue of 
biodiversity measured by biodiversity FP. 
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Nomenclature 

Superscripts 

LI harvesting and supply layer 

L2 collection and pre-processing layer 

L3 main processing layer 

L4 use layer 

ind indirect 

road road conditions 

sub substituted product 

tr transport 

Sets 

F set of footprints 

/ set of supply zones 

J set of demand locations 

M set of collection and intermediate process centres 

N set of process plants 

P set of products 

PS set of substituted products 

T set of technology options 

Subsets 

PD set of directly used products (subset of P) 

PI Set of intermediate products (subset of P) 

PIP = PI x PP set of pairs of intermediate product and produced 
product (if a produced product is produced from 
a given intermediate product) 

PP set of produced products from plants (subset of P) 

PSP = PSxPP set of pairs of substituted product and produced 
product (if a produced product substitutes 
conventional product) 

PT=PIxT set of pairs of intermediate product and applicable 
process technology for it 

Indexes 

f index for footprints 

i index for supply zones 

j index for demand locations 

m index for collection and intermediate process centres 
n index for process plants 

p index for products 


pd index for directly used products 

pi index of intermediate products 

pp index for directly used products 

ps index for substituted products 

t index for technology options 

Scalars 

A supply-chain network’s total area, km 2 

w H2 o mass fraction of water 

Parameters 

ei'Aj specific environmental footprint/of intermediate 

product pi at layer a,a= {1,2}, t/t or GJ/t or km 2 /t 
eij+r/ specific environmental footprint/of intermediate 

product pi and the selected technology t at processing 
layer, t/t or GJ/t or km 2 /t 

eij;^ specific environmental footprint/of product p at use 
layer, t/t 

eip“^ specific environmental footprint/of substituted product 
ps, (t or GJ)/(MWh or GJ or t) 

/p° nv ,L2 conversion factor of intermediate product pi by pre¬ 
processing 

Jp“pp substitution factor between use of conventional energy ps 

and biomass energy pp 

Continuous variables 

ENVB}“ direct environmental footprint/at layer a, t/(km 2 y) or 
GJ/(km 2 y) or km 2 /(km 2 y) 

q' n pi L ’ production rate of intermediate product pi at supply zone 

TT7I3 ' ,t/y 

Q’n pi t ' mass-flow of intermediate product pi to the selected 

technology t at the process plant n, t/y 
q x!y L p’ Lb mass-flow of product p from object x in layer a to object y 
in layer b, t/y 

D*|y Lb distance between object x in layer a and object y in layer 
b, km 

goad,La,Lb roac ] conc ]jtion f ac t 0 r between object x in layer a and 
object y in layer b 

eipj a,Lb transport environmental footprint/from layer a to the 
layer b, t/(t km) or GJ/(t km) 

Inpippt 3 mass-flow of produced products pp from intermediate 
product pi with the selected technology t at the process 
plant n, t/y 

ENVB| nd indirect environmental footprint type/, t/(km 2 y) or 
GJ/(km 2 y) 

ENVBf total environmental footprint type/, t/(km 2 y) or 
GJ/(km 2 y) or km 2 /(km 2 y) 

Abbreviations 

CC Collection and pre-treatment centre 

CFP Carbon footprint 

DDGS Distiller’s dried grains with solubles 

EFP Energy footprint 

FEFP Food-to-energy footprint 

FP Footprint 

GHG Greenhouse gas 

LCA Life cycle assessment, also Life cycle analysis 

LFP Agricultural land footprint 

MCO Multi-criteria optimisation 

MINLP Mixed-integer non-linear programming 

MOO Multi-objective optimisation 

MSW Municipal solid waste 

RFP Relative footprint 

WFP Water footprint 

WPFP Water pollution footprint 
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